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In this and subsequent papers, one dimensional system of Dirac fermions with a random- varying 
mass is studied by the transfer-matrix methods which we developed recently. We investigate the 
effects of nonlocal correlation of the spatial-varying Dirac mass on the delocalization transition. 
Especially we numerically calculate both the "typical" and "mean" localization lengths as a function 
of energy and the correlation length of the random mass. To this end we introduce an imaginary 
vector potential as suggested by Hatano and Nelson and solve the eigenvalue problem. Numerical 
calculations are in good agreement with the results of the analytical calculations. We obtain the 
relation between the localization length of states and the correlation length of the random mass. In 
subsequent paper, we shall study Dirac fermions with a random mass of long-ranged correlations. 



PACS: 72.15.Rn, 73.20.Jc, 72.10.Bg 
I. INTRODUCTION 

One of the most important problems in condensed 
matter physics is the localization phenomena in random- 
disordered systems At present it is believed that al- 
most all states tend to localize in disorder systems in two 
and lower dimensions. In some special cases, however, 
some specific states remain extended even in the pres- 
ence of strong disorders. System of Dirac fermions with a 
random- varying mass(RMDF) in one dimension has been 
studied from this point of view §-|l. RMDF IS a con- 
tinuum field theory of the random hopping tight binding 
model and it is related with the random XY spin chain, 
the random Ising model, polyacetylene system, etc. 

In most of studies on the random-disordered systems, 
only the short-range white-noise random variables are 
considered. In the previous papers |lC|-|l^ we studied 
the effect of nonlocal correlation of the random mass 
on the extended states which exist near the band cen- 
ter. For numerical studies, we reformulate the system 
by transfer-matrix formalism, and obtained eigenvalues 
and wave functions for various configurations of random 
telegraphic mass 0. We verified that the density of 
states obtained by the transfer-matrix methods is in good 
agreement with the analytical calculation by supersym- 
metric methods in Ref. [ pd] |. In this paper we shall in- 
troduce an imaginary vector potential into the system of 
the RMDF and study a localization-delocalization phase 
transition by varying the magnitude of the vector poten- 
tial. Through this study we shall obtain the localization 
length of the states quite accurately as a function of en- 
ergy and the correlation length of the random mass. This 
method of calculating the localization length is based on 
the idea by Hatano and Nelson and it is totally dif- 
ferent from the methods which have been used so far [f4| . 



In this and subsequent papers we shall study effects of 
long-range spatial correlations of the random mass, es- 
pecially on dclocalized states, localization lengths, etc. 
This work is partially motivated by the recent studies on 
the one-dimensional Anderson model with random po- 
tentials with long-range correlation. There it is shown 
that there exists a mobility edge at a finite energy if the 
correlation is long-range enough. As we explained above, 
extended states exist at the band center in the present 
model even for the white-noise random mass. Then it 
is interesting to see if additional extended states appear 
as a result of the long-range correlation of the random 
mass. 

This paper is organized as follows. In Sec. 2, we first 
explain the imaginary- vector potential methods (IVPM) 
for calculating the typical and mean localization lengths. 
Then we compare numerical results with the analytic cal- 
culations in order to verify the validity of the IVPM. 
In Sec. 3, we shall consider the system of the RMDF in 
which a Gaussian distribution is employed for distances 
between (anti-)kinks or jumps of the telegraphic random 
mass. We calculate the correlation length of the random 
mass and also the localization lengths as a function of the 
mean distance and standard deviation. Then we obtain 
relation between the correlation length of the random 
mass and the localization lengths. Section 4 is devoted 
to conclusion. 

The studies on the Anderson model etc., in one di- 
mension employ specific techniques like the Hamiltonian 
mapping, the renormalization group, etc. Compared 
with them, our method are rather straightforward. Be- 
fore calculating the localization lengths, we obtain land- 
scapes of the wave functions which are useful to get phys- 
ical pictures of the delocalization transition. 
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II. MODEL, TMM AND IVPM 

We shall consider a Dirac fermion in one spatial di- 
mension with a coordinate-dependent mass to (a;) and in 
an imaginary vector potential g, whose Hamiltonian is 
given by, 



n 

h 



ia^'idx +g) + m{x)(T'^ , 



(1) 
(2) 



where a are the Pauli matrices. This fermion model is a 
low-energy effective model of random- hopping tight bind- 
ing models, random-bond spin chains, etc |^]. We intro- 
duce the components of as — {u,v). In terms of 
them the Dirac equation is given as. 



d_ 

dx 
d_ 

dx 

(We follow the notations in Ref. Q.) From Eqs. 
obtain the Schrodinger equations. 



+ g + m{x) j u{x) — Ev{x), 
— g + m{x) \v{x)—Eu{x). 



(3) 



dx'^ 



2.9 



dx 



m'{x) + {rn?{x) — 9^)j '"(^) 
= E'^u{x), 



(4) 



and similarly for v{x). 

In this paper we restrict the shapes of m[x) to multi- 
kink-anti-kink configurations [l5| . The multi-soliton- 
antisoliton configurations are given by. 



•m{x) 



m{9{x — ai) — 1) 



^to(0(-x + /3,)-1), 



(5) 



where ai's(/3j's) are positions of kinks(anti-kinks) and 
they are random variables in the present system. An 
example of to(x) is given in Fig.l. 



' ■m(x) 



FIG. 1. An example of configurations of solitons and 
anti-solitons. 



If we vary the distances I between soliton and anti- 
soliton according to the exponential distribution like. 



(6) 



where A is a parameter, then m{x) has the following cor- 
relation ||], 

A 

[ m{x) m{y) ]e„s = exp {-\x-y\/X), (7) 
A 



where 



essentially corresponds to the height of the 



soliton and anti-soliton, i.e., to in (||). From (^, A is 
the correlation length of the random mass and the limit 
A ^ corresponds to the white-noise case. In the subse- 
quent papers p^ , p7[ | , we shall study Dirac fermions with 
long-range correlated random mass by using the meth- 
ods examined in this paper. There we expect some in- 
teresting phenomena like existence of nontrivial mobility 
edge, nonuniversality of the multi-fractal exponents, etc. 
Studies in this section show that the IVPM for calcu- 
lating the localization lengths are reliable and we shall 
use them for studies on random systems with long-range 
correlated disorder. 

For the vanishing imaginary vector potential, we can 
solve the Schrodinger equations in (^) under the pe- 
riodic boundary condition with various multi-soliton- 
antisoliton configurations of 'm{x) using the transfer- 
matrix method(TMM). We can also obtain the energy 
spectrum and wave functions ]To| using this method. 

We review the TMM method briefly. Let us consider 
the Dirac equation with a mass configuration which 
has a kink at a; = 0, 



m{x) = m{29{—x) — 1) -I- toq, 
where 9{x) is the Heaviside function, 



and it is assumed that 



(x < 0) 
{x > 0) 



(8) 



(9) 



TO > Too > 0. 

Substituting (|) into @, we obtain a Schrodinger equa- 
tion with a potential which is a combination of the delta 
function and the step function, 

V{x) — ^m'{x) + m'^{x) — ±2m(5(x) -I- 4tooto0(— a;) 



+ (to — TOo)^ 



(10) 



For such a field equation with a delta-function-type po- 
tential, the continuity of the wavefunction u{x) leads au- 
tomatically to the conditions on its derivative u'{x). For 
example, if we rewrite the Schrodinger equation as 
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+ 2m S{x) + Amom 6{—x)j u{x) 

^ {E^ ~ {m-mof) u{x), (11) 



(here we restrict g to zero) and integrate Eq.(^l]) with 
respect to x in the range — e < x < e, we have 

.0 

-u'(0 + e) + u'{0 - e) + 2fh u(0) + -imofh / dx u{x) 
= {E^ - (to - mof) [ dx u{x). (12) 



We see the r.h.s. of Eq.(12) vanishes in the hmit e ^ 
due to the continuity of u{x). A similar argument appHes 
to v{x) and we have conditions, 



u'(0 + 0) - u'{0 - 0) = 2mu(0), 
v{0 + 0) - v'{0 - 0) = -2tow(0). 

We look for solutions under these conditions. 
Let the energy eigenvalue E satisfy 

m — Too > -E > 0, 

then u{x) takes the form, 



u{x) 



A 



B e 



C e'"^"' + D e 



(x <0), 
(a;>0). 



Ki = \/{m + TOo)2 - E'^, 

K2 = \/(to - TOo)2 - E'^. 



From the continuity of m(x), 

A + B = C + D, 

and Eq.(p^) gives 

K2(C -D)~ Ki{A - B) = 2m{A + B). 

Hence, 

u{x) = Ae"^"" + Se^"!^, for x < 
and for a; > 0, 

2to + Ki + K2 2to — Ki + K2 



(13) 



(14) 



(15) 

(16) 
(17) 



(18) 
(19) 
(20) 



u{x) — 



2k2 

2m — Kl + K2 



-A 



■A- 



2k2 

-2to + Ki + K2 



-B 



B 



(21) 



2k2 2k2 

If we let Too = 0, ki = K2 = k, Eq.(pl|) simplifies to 



m(x) = < 



m + K , TO „ 
A H B 



TO 



-A 



-m + K 



B 



(x < 0) 
(x > 0) 

(22) 



Values of A, B and E are determined by boundary con- 
ditions, e.g., u{~L/2) = u{L/2) = for a given system 
size L. 

We restrict toq to zero in the rest of the discussion. 
Since the wavefunction u{x) is expressed everywhere as 



u{x) = A e'^^ + B e" 



(23) 



we can represent the eigenfunction in terms of coefficients 
A and B. By using matrix representation and from ([2l]), 
the conditions ( p^ ) give the following relation between 
the coefficients A and _B, 



A{x > 0) 
B{x > 0) 



T 



A{x < 0) 



B{x < 0) 
/ 1 + — — \ 



T 



V 



(24) 



TO TO 
— 1--J 



For kink instead of anti-kink, we should replace to with 
(— m) in Eq.(p^) and we obtain 



A{x > 0) 
bIx > 0) 



T- 



A{x < 0) 
B{x < 0) 



(25) 



Let us define 



R{k, a) 







m 



(26) 



and TM for the configuration of an anti-kink and a kink 
is given by 



R{b)T-^R{a)T 

— r?th\ ( ~ "^^"^ sinhfto 2(j){l — (j)) sinh na 

" ^ ' \ 20(1 + (t>) sinh Ko e-""" + 20^ sinh Ka 

= T(a,6), (27) 



where a is the distance between the kink and anti-kink. 

We impose periodic boundary condition on the wave- 
function. The above argument can be generalized to an 
arbitrary number of kink and anti-kink pairs readily and 
we have the following equation expressed by the transfer 
matrices. 



det 



T{eJ)---T{c,d)T{a,b)- 



1 
1 



= 0. (28) 



We are not able to solve ( p8| ) analytically for arbitrary 
configuration of pairs of kinks. However it is not so dif- 
ficult to solve it numerically. We can also easily obtain 
eigenfunctions after having energy eigenvalues by using 
TM method. It is straightforward to extend the above 
formalism for nonvanishing toq. 
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Let us introduce the constant imaginary vector poten- 
tial g. Effect of the imaginary vector potential was dis- 
cussed by Hatano and Nelson Let us denote the 
eigenfunction of energy E for g = as ^Jq{x), and sup- 
pose the shape of ^0(2;) as 



*o(a;) = exp 



(29) 



where ^0 is the localization length and Xc is the center of 
this localized state. When we turn on the constant imag- 
inary vector potential g, the eigenfunction is obtained 
from ^'0(2:) by the "imaginary" gauge transformation, 



= exp 



^0 



g (x- Xc) 



(30) 



This means that the localization length of this eigenfunc- 
tion is 

Co 



^ ~ 9t,o 



(31) 



At the point g = the localization length for x < Xc 
diverges. If the imaginary vector potential g gets larger 
than this value, for x < x,, becomes negative, and the 
eigenfunction cannot exist at the same energy since it 
cannot satisfy periodic boundary condition (Fig. 2). 

So if localized eigenstate with energy E disappears at 
g = gc as g is increased, then the localization length of 
the eigenstate 'ifo(x) is 1/gc- Actually it is shown that for 
g > gc, energy eigenvalue of the state has an imaginary 
part and the state is extended as we shall see shortly . 

The TMM can be easily extended for the case of nonva- 
nishing g. We obtain energy eigenvalues and eigenfunc- 
tions numerically. In actual calculation the system size 
is finite, and the energy eigenvalues of the wavefunctions 
for 5 = and 5 7^ are slightly different as we mentioned 
above. We must trace the changes of energy eigenvalues 
of states for nonzero g on the two dimensional {E — g) 
plane and find the end point of the locus (Fig. 3). 

In the numerical calculation, we first find the energy 
eigenstate for g — 0. We denote its energy eigenvalue 
as -Bo. Then we increase g and search the eigenstate in 
the region Eq - SE < E < Eq + SE. (We set SE very 
small.) Let us suppose that the eigenstate vanishes in this 
energy region at the value g = gi- We denote the energy 
eigenvalue of the eigenstate for g — gi as Eg^ . We change 
the energy region for searching eigenstate as Eg^ — SE < 
E < Eg^ + 6E and let g larger than gi . We repeat this 
procedure. If we cannot find energy eigenvalue at the 
value g larger than gc, we judge that the localization 
length of the eigenstate for g = is l/gc- 



g=o 






X 



'P(x) X 



This wavefunclion cannot 
satisfy periodic boundary 
condition. 



FIG. 2. The imaginary vector potential method 
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J Search range (2nd step) 


Egi-5E| \ 
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Eo-5E 


' \Eo+5E Energy 




^Eo 



FIG. 3. The procedure of imaginary vector potential 
method in the numerical calculation 

First of all, we show that the delocalization transition 
actually occurs at finite g. In Fig. 4, we show the wave 
functions of a low-lying state in vanishing and nonvan- 
ishing imaginary vector potential. For g = 0, the state 
is obviously localized whereas at g = 0.03 the state be- 
comes extended and the energy eigenvalue has an imag- 
inary part. As discussed in Ref. the density distri- 
bution of a particle is given by {x , ~ g)'^ {x , g)\ where 
'i'{x, —g) is equal to the left eigenfunction. 
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l1'-g(x)'Pg(x)l 



If-gCxjfgCx)! 



E=0.1014+0.00657i 



FIG. 4. An example of the localized and the extended wave 
functions in vanishing and nonvanishing imaginary vector po- 
tential 3 = and g = 0.03. The state of jr = 0.03 has a 
complex energy eigenvalue. 



Let us turn to the localization length. For the white- 
noise case [m(x)m{y)]ens = ^ S{x — y), "typical" local- 
ization length or the inverse of the Lyapunov exponent 
was obtained as P,n8,n9[ 



^t{E) = \\nE/2A\. 



(32) 



Numerically the typical localization length £,t{E) is ob- 
tained by averaging over localization lengths of all eigen- 
states with energy E. On the other hand, Balents and 
Fisher calculated the averaged Green function and ob- 
tained the mean localization length from the spatial de- 
cay of the Green function jsj . The result is 



U{E) = \\nE/2A\' 



(33) 



Case of nonlocally correlated random mass was studied 
in Refs. ||ll| , |l^ and S,m{E) is obtained as a function of A 
in Eq.(0). 

By numerical calculation we obtain both the typical 
and mean localization lengths. As we mentioned above, 
the typical localization length is the average over all solu- 
tions of the Schrodinger equation (^), whereas the mean 
localization length is determined by the states which 
make dominant contributions to the Green function, i.e., 
which have large localization length. 

The result of the numerical calculation of the typical lo- 
calization length of the white- noise case is given in Fig. 5. 
We show the ratio of the numerical results to the analyt- 
ical expression in Eq. ( ^2| ) in order to compare these two 
results. Therefore if the energy dependences of the lo- 
calization lengths obtained numerically and analytically 
are the same, this ratio should be constant. In Fig. 5, the 
ratio seems constant over the whole range of E. 

In Fig. 6, we show the numerical results of the "mean" 
localization length. Here we use the solutions to the 



Schrodinger equation which have long localization length. 
More precisely, "large" localization length ^ means 
the one which satisfies ^ > (the "typical" localization 
length) -|- 1.0 a in each energy slice. From Fig. 6, we con- 
clude that the energy dependence of the mean localiza- 
tion length obtained numerically is in agreement with 
Eq.(|3|). 

Typical Localization Length 




0.05 0.1 0.15 0.2 0.25 

Energy 

FIG. 5. The comparison between the analytical and the 
numerical results of the typical localization length for the (al- 
most) white-noise case, i.e., small- A case: In the numerical 
calculation, we set L(system size)= 50, A = 1/2, A — 1/60 in 
Eq.(0), and the energy slice AE — 0.01. The ratio is normal- 
ized at i5 = 0.15, and this result is averaged over about 15000 
eigenstates. The analytical and the numerical results are in 
good agreement. 



Mean Localization Length 



Mean Localization Length / lln E/4AI 



i'.JlH!l.li..«i.i;„...y.» 



Mean Localization Length / lln E/4AI 



0.05 0.1 0.15 0.2 0.25 

Energy 

FIG. 6. The comparison between the analytical and the 
numerical results of the mean localization length: Here we 
picked up the eigenstates which have large localization length. 
(See the text.) We show the ratios of the numerical calculation 
both to Eqs.dl) and (||). We set L = 50, A = 1/2 and 
the energy slice AE — 0.02. The ratios are normalized at 
E = 0.15. 



We shall turn to the case of the nonlocally-correlated 
disorder. In the white-noise limit (A = case in Eq.(^), 
the localization lengths diverge only at i? = 0, that is, 
extended states exist only at E = 0. If we let A > 0, 
the random mass becomes nonlocally-correlated, and the 
critical energy or the mobility edge at which the delocal- 
ization transition occurs may change. 

We investigate the "typical" and "mean" localization 
lengths in the case of nonvanishing A's. The behaviour of 
the "typical" and "mean" localization lengths obtained as 
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in the white-noise case are given in Figs. 7 and 8. It seems 
that there is no A-dependence in the typical locahzation 
length. On the other hand, Fig. 8 shows that the mean 
localization length has a small but finite dependence on 
A. 

From the above calculations, we conclude that the ef- 
fect of the short-range correlations in disorders is not so 
large. Especially the result indicates that the delocaliza- 
tion transition occurs at £^ < 0.03. (If the mobility edge 
exists at Ec > 0, the ratio in Fig. 7 or 8 must diverge 
at Ec-) The delocalization transition probably occurs at 
E^O. 

Typical Localization Length 

2.5 I 1 1 ^ 1 1 

J, = 1/60 ' ' — ' 

1=1/4 



1.5 ■ 

(2 1 E* * ■ g ' tw i ll 1 1 ti i * *i S * 

0.5 ■ Typical Localization Length / lln E/4AI . 

I ' • • • 

0.05 0.1 0.15 0.2 0.25 

Energy 

FIG. 7. The behaviour of the localization length in the case 
of nonvanishing A's: Here we show the ratios of the typical 
locahzation lengths to |ln£;/4yl| (Eq.(||)). We set L = 50, 
A = 1/2 and the energy slice AE = 0.01. The ratio is nor- 
malized a.t E = 0.15. 



Mean Localization Length 
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1=1/4 
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Mean Localization Lenj 


»th / lln E/4AI- 



0.05 ai 0.15 0.2 0.25 

Energy 

FIG. 8. The behaviour of the localization length in the case 
of nonvanishing A's: We show the ratios of the mean locahza- 
tion lengths to | \nE/AA\^ (Eq.(||)). We set L = 50, yl = 1/2 
and the energy slice AS = 0.02. The ratio is normalized at 
E = 0.15. 

In the previous paper we calculated the localiza- 
tion length for the random mass with the short-range cor- 
relation (0) . We obtained the "mean" localization length 
to the 1st order of A by means of the Green function 
method. The result is 

m = ^(i^ + ^A^i^) +0(P). (34) 

In Fig. 9 we show the ratios of the numerical result to the 
analytical calculation up to the 0th and the 1st order of A. 



This shows that the analytical result with the 1st order 
correction of A is in better agreement with the numerical 
result, but the correction is small. 



Mean Localization Length 



1 St order 












0th order 





0.05 0.1 0.15 0.2 

Energy 

FIG. 9. The comparison between the numerical result and 
the analytical one in the case of relatively large A: We show 
the ratios of the mean localization length calculated numer- 
ically to the ones obtained analytically (Eq.(p^)). Here we 
used the analytical results of the 0th and the 1st order of A. 
We set A = 1/4, A = 1/2 and L — 50. The energy slice is 
AE — 0.02, and the ratio is normalized at i5 = 0.15. 

From the investigations given as far we can conclude 
that the numerical methods used in this paper are reli- 
able for calculating the localization lengths. It is very 
interesting to study the case of disorders with a long- 
range correlation. We expect that a nontrivial mobility 
edge Ec > exists for a certain long-range correlated 
random mass. Actually the one-dimensional Anderson 
model with long-range correlated disorder was studied 
, and it is shown that there exists a nontrivial mobil- 
ity edge. 

We can also calculate exponents of the multi-fractal 
scaling by the TMM [|l^ . These problems are under 
study and results will be reported in a future publication 




III. GAUSSIAN-DISTRIBUTED DISTANCES 
BETWEEN KINKS: CORRELATION LENGTH 
VS. LOCALIZATION LENGTHS 

In subsequent papers ||l^,|l3 we shall study the local- 
ization length with long-range correlated disorder with 
the TMM and IVPM. Before that study, we shaU consider 
the case of Gaussian distribution for distances between 
kinks in this section. That is, the probability distribu- 
tion function Pq for distance I is given by the Gaussian 
distribution, 

Pg{1) - cxp{-(/ - A^)V2an, (35) 

where fi is the mean distance and a is the standard devi- 
ation. If a is small, eigenstates at low energy (eigenstates 
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which have a small number of nodes) spread over whole 
system, and a mobility edge may exist at a nonzero en- 
ergy. Physically, the mean distance corresponds to the 
kink density ( it can be viewed as impurity density in the 
related models such as the random hopping tight binding 
model, spin chain, etc.) and standard deviation controls 
the degree of randomness. This is in sharp contrast to 
P{1) in Eq.(|^) in which only one parameter A exists. 

We first show the relation between localization lengths 
of eigenstates and the deviation a. The result is shown 
in Fig. 10. There are several peaks of localization length 
in these figures with small a. We know that the energy 
spectrum is discrete and all states become extended in 
the case of periodic m(a;), i.e., the limit ct — > 0, as dic- 
tated by the Bloch theorem. For small but finite tr, we 
expect that energy spectrum is not discrete and eigenval- 
ues spread within narrow energy band. We also expect 
that localization lengths of states are very large, but fi- 
nite for nonvanishing a. 



fL=l , (1=0,050, L=50 



H=l , 1=0.250. L=50 
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FIG. 10. Localization length and standard deviation of 
Gaussian distribution: Here we set L(system size)= 50. In 
the case of small a, the points concentrate around the eigen- 
values of periodic system and they form several peaks. As 
we let cr larger, the peaks disappear and these figures become 
similar to the one in case of white-noise disorder. (Extended 
state exists only at E = and other states localize.) 



In a system of finite size, we regard state as extended 
one if the localization length is larger than the system 
size. We can see that localization lengths of several states 
become larger than the system size if we take a small. 
This means that some eigenstates become extended in 
finite size system. 

If we let a much larger, these peaks of localization 
length get smaller and disappear except for the one at 
E ^ 0. We recall that the state at i? = is "isolated" 
extended state in the system with white-noise correlated 
random mass. In the case of large cr, the random mass 
m{x) becomes short-ranged and therefore all states ex- 
cept those close to E' = tend to localize. 

Next we show the two-point correlation function of 
random mass [m(x)m(0)]ens (Fig-H)- It oscillates due 
to Gaussian distributed randomness. The period of this 
oscillation reflects mean value of Gaussian distribution. 
The value of envelope of [m(a;)m(0)]ons approaches to 
zero for larger x, and the envelope shows exponential de- 
cay, as we can see in Fig.ll. This behaviour is the same 
as in the case of exponentially distributed kink distances. 
In the present case, however, we can obtain the system 
with long-range correlated disorder in which there are a 
large number of kinks. (In the case of exponential distri- 
bution P{1) in Eq.(^, there are a small number of kinks 
in the system with large correlation length, and the distri- 
bution with P{1) is not sufficient for studies on practical 
systems like random spin chains.) The correlation length 
of random mass is very large and exceeds system size for 
small cr. 



1 

o.s 

0.6 
0.4 
0.2 


-0.2 
-0.4 
-0,6 
-0,8 



Correlation of m(x) 



0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0,45 0,50 
a (/system size) 
Correlation of m{x) 





-0,5 



0.05 0,10 0.15 0.20 0,25 0.30 0,35 0.40 0,45 0,50 0,55 
a {/system size) 



FIG. 11. Correlation of random mass with Gaussian dis- 
tributed randomness: Here we set ^ = 1, a = 0.20 and 
//(system size)= 50. We show the behaviour of envelope in 
the logscale plot. From this figure, we know that m{x) is 
exponentially correlated in this case. 
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The relation between a and the correlation length A 
(decay rate) of random mass is shown in Figs. 12 and 13. 
(A is defined from the envelope of correlation of m(x) 
as [m(a;)m(a; + a)]ens ~ exp(— a/A).) From Fig. 13, we 
conclude that these quantities are related as 

A = Aa-^. (36) 

In the case shown in Fig. 13, A and v are estimated as 



A = 4.91 X 10"^, V = 1.881 {L = 50, /i = 1), 
A = 2.67 X 10"^, V = 1.837 {L = 100, [i = 1), 

where L is the system size. 

a and Correlation Length of m(x) 



(37) 



a 1 
e 

% 0.8 

I 0.6 

I 

I 



L=50, (1=1 
L=100, |1=1 



0.05 0.1 0.15 0.2 0.25 03 0.35 0.4 

a 

FIG. 12. Relation between deviation a of Gaussian dis- 
tribution and correlation length A of random mass(l): A is 
defined as [m{x)m{x + a)]cns ~ exp(— a/A). In the case of 
small (T, correlation length of m{x) is much larger than the 
system size. 



o and Correlation Length of m(x) 



E 

~ -1 



L=50, n=l 

(Y=-5.3I5-1.S8IX) 

L=100, |i=l 

(Y=o.922-1.837X) 



S a (Logscale) 

FIG. 13. Relation between deviation a of Gaussian distri- 
bution and correlation length A of random mass(2) (Logscale 
vs. logscale plot): Points for fixed system size are almost on 
one line. 

We observe the typical localization length of the states 
at a fixed energy and also its dependence on system 
size. We have already seen that energy eigenvalues of 
the states concentrate around several energy values. In 
order to calculate typical localization length, we pick up 
the states around second peak from the lowest energy and 
average localization lengths of these states. (See Fig. 10. 
In the actual calculation, we picked up the states around 
E ~ 0.1 in the system with L = 50, ~ 1, and the 



ones around half(quarter) of this energy value in the dou- 
ble(four times) size system.) 

The dependences of typical localization length on a 
and system size are shown in Fig. 14. We expect that 
the a dependence of localization length ^ can be fit by 
following function. 



e = A |ct - CTc 



(38) 



where CTc is critical value and rj is critical exponent. 
Strictly speaking, this relation holds true only in infi- 
nite size system. We expect that we can assume this 
parameterization if the system size is large enough. 

The result of fitting is shown in the Table 1. We con- 
clude that the choice of this fitting function is valid from 
these values and also from the fact that the critical 
exponents for different system sizes are almost the same. 



(J dependence of Typical Localization Length 




U 3-5 
3 

1,2.5 
c 

2 2 
o 

S 1.5 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



o dependence of Typical Localization Length 



L=100 (2a error) 

3 

I 2.5 

J 2 
a 

I 1-5 

I ' 
^ 0-5 

"I 

K 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



1^ a dependence of Typical Localization Length 



|>2.5 








0.05 0.1 0.15 0.2 0.25 0.3 0.35 04 



FIG. 14. Typical localization length and deviation a: We 
see the a dependences of typical localization lengths for sec- 
ond peak in various size systems. Here we set /i(mean 
distance) = 1. We also show the results of fitting. (They 
are shown as the curve in each figure.) 
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Table 1 

L AjxlO--) a.(xlO-^) , 

50 12.03 ±3.41 4.57 ±10.98 1.416 ± 0.057 0.493 
100 7.62 ± 1.83 3.70 ±8.88 1.411 ±0.045 0.479 
200 4.09 ±0.76 1.23 ±6.59 1.386 ±0.031 0.754 



The coefficient A becomes smaller as we make the sys- 
tem size larger. This means that the states become more 
localized in larger system, cjc also approaches zero as 
the system size is going to infinity. Energy eigenvalues 
of states around second peak approach zero in larger sys- 
tem. From these facts, we expect that all states at i? > 
localize and there exist isolated extended states only at 
E = in infinite size system with nonzero a. This is the 
same situation as in the case of white-noise disorder. 

From the above calculations, we can obtain the relation 
between the correlation length of the random mass A and 
the localization length ^ as follows, 

e-CA^ (39) 

from Eqs.(^6|) and (^). The value of r is estimated as 
around 0.8. 

In the random bond XY model, the localization length 
^ is related with the correlation length of spins. On 
the other hand, the impurity density and/or aperiodicity 
of the positions of impurities control the parameter A. 
Therefore we hope that the result ( |39| ) can be observed 
at least qualitatively in the random spin chains. 

IV. CONCLUSION 

In this paper we studied the random-mass Dirac 
fermion in one dimension by using the TMM and IVPM. 
In the first half, we explained and verified the validity 
of the IVPM by comparing the obtained results with the 
analytical calculations. In the second half, we investi- 
gated the effects of the nonlocal correlations of the ran- 
dom mass by using Gaussian distribution for distances 
between jumps of m{x). We obtained relation between 
the correlation length and localization length. In the sub- 
sequent papers [|l^,0 , we shall study the Dirac fermion 
with long-range correlated random mass. 
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